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Some recent results on evaluating Feynman integrals 

V.A. Smirnov''* 

Nuclear Physics Institute of Moscow State University, 
Moscow 119992, Russia 

Some recent results on evaluating Feynman integrals are reviewed. The status of the method based on Mellin- 
Barnes representation as a powerful tool to evaluate individual Feynman integrals is characterized. A new method 
based on Grobner bases to solve integration by parts relations in an automatic way is described. 
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1. Introduction 

Perturbativc quantum-theoretical amplitudes 
are expressed in terms of Feynman integrals over 
loop momenta. Usually, one turns immediately 
to scalar Feynman integrals using some projec- 
tors and obtains a family of scalar Feynman in- 
tegrals with the same structure of the integrand 
and various powers of propagators (indices): 
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where fc^, i = l,...,/i, are loop momenta and 
the denominators Er are either quadratic of lin- 
ear with respect to ki and external momenta 
qi, . . . , QN ■ By default, the integrals are dimen- 
sionally regularized with d = 4 — 2e. 

If the number of Feynman integrals needed for 
a given calculation is small or/and they are sim- 
ple, one evaluates, by some methods, every scalar 
Feynman integral of the given family. Various 
methods are used, in particular, alpha and/or 
Feynman parameters, Mellin-Barnes (MB) rep- 
resentation |1I2| and differential equations In 
the next section, the method of MB representa- 
tion is briefly reviewed. 

If it is necessary to evaluate a lot of compli- 
cated Feynman integrals ^ the standard way is 
to apply integration by parts (IBP) ^ relations 
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in order to construct an algorithm that gives 
and expression of a given Feynman integral as 
a linear combination of some master integrals, 
F(ai,...,a„) = Y.iC{ai,- ■ ■ ,a„)Ii Here pi = 
ki,...,ph = kh,Ph+i = gi, • ■ • ,Ph+N = qN- Once 
one has a solution of this (reduction) problem, 
it becomes sufficient to evaluate only master in- 
tegrals. The IBP relations can be written as 
fi ■ F where fi are polynomials in the shift op- 
erators Yr^^ defined by {Y^^^ • F)(ai, . . . , a„) 
F{ai, . . . ,ai + 1, . . . , an). 

The first attempt to make the reduction pro- 
cedure systematic was based on the fact that the 
total number of IBP equations grows faster than 
the number of Feynman integrals satisfying the 
condition \ai\ < M at M ^ cxo so that, at suf- 
ficiently large M , one obtains an overdetermined 
system of equations |5I6| which can be solved. 
There is already a public implementation of this 
algorithm on a computer |2] . Another attempt |H] 
is based on a special parametric representation. 

Other attempts are based on the use of Grobner 
bases (GB) This idea was first suggested in 
|10j . where IBP relations were reduced to differ- 
ential equations. To do this, it is assumed that 
there is a non-zero mass for each line. The typical 
combination aiYi, is then naturally transformed 
into the operator of differentiation in the corre- 
sponding mass. An attempt to use GB associated 
with the shift operators was made in [llj . How- 
ever, the corresponding algorithms can now work 
only in simplest cases, for n = 2. 

In Section 3, another approach ^21 based on 
GB is briefly characterized. In conclusion, I shall 
discuss some perspectives. 
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2. Evaluating by MB representation 

The method of MB representation was success- 
fully applied to evaluate massless on-shell double 
|1I2I13I14I151 and triple |16I17| boxes, with results 
written in terms of harmonic polylogarithms 18j , 
double boxes with one leg off shell ^6.19^ and mas- 
sive on-shell double boxes |2()I21| The method is 
based on the MB representation 

1 _ r+'°^' r(A + z)r(-z) dz 

applied to replace a sum of two terms raised to 
some power by their products in some powers. 

Experience shows that a minimal number of 
MB integrations is achieved if one introduces MB 
integrations loop by loop, i.e. derives a MB repre- 
sentation for a one-loop subintegral, inserts it into 
a higher two-loop integral, etc. For example, for 
the tennis court graph shown in Fig. 1, one can 
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start from a MB representation for a lower box 
subintegral with three legs off shell, then insert it 
into the double-box subintegral to obtain a MB 
representation for the double box with two (up- 
per) legs off shell, insert it into the whole graph 
and obtain an eightfold MB representation for the 
given family of integrals with general powers of 
the propagators and the power of the numerator 
[(^1 + : where /1.3 are the momenta flow- 

ing through lines 1 and 3 in the same direction: 

T(ai,...,aii;s,t;e) 

(i7r'*/2)'(-lf (-s)8-«-3^ 

0^=2,4,5,6,7,8 r(aj)r(4 - 04567 " Se)^^ 



r(a - 8 -I- 3e -I- w)T{S - a - ie - w) 
r(ai - Z2)r(a3 - 2;3)r(a9 - Z6)r(aio - 247) 
r(a5 + Zi4)r(a2 -I- 2:56)r(2 - w + Z5) 
^ r(4 - ai23 - 2e + zi23)r(8 - a - - z^) 
r(2 - 0457 - £ - 2i3)r(aio -2 + w - Z457) 

r(ai234567,ll - 4 -I- 2e -I- 2:4567) 
xr(2 - 0567 - e - 2;i24)r(a4567 - 2 + e Z1234) 

xr(2 - 023 - e + 2:13 - 25)r(a9 -2 + w - z^q) 
xr(2 - ai2 - e + Z12 - Z567)r(a7 + Z123) 
Xr(ai23 - 2 + e - Z123 + Z567)r(z57 - Zi) 
xr(4 - 089,10 - e - w + Z4567) , 

with ai2 = ai + a2,a = ai...^ii, Z14 = zi-\- Z4, etc. 

Deriving MB representations for general indices 
is very useful because one can apply it for various 
partial cases and obtain crucial checks. For ex- 
ample, this very representation was used in |17) to 
calculate r(l,...,l,— 1), ina Laurent expansion 
in e and to check cross order relations in = 4 
SUSY gauge theories [22]. Let us stress that one 
can check such a cumbersome representation in an 
easy way by considering two partial cases: when 
one contracts horizontal lines, i.e. in the limit 
01,03,04,06,08 0, or vertical lines, i.e. at 
02,05,07,09,010 — ^ 0. In both cases, one obtains 
recursively one-loop integrals which can be eval- 
uated in terms of gamma functions for general 
e. On the other hand, taking such limits reduces 
to calculating residues in some integration vari- 
ables. Consider, for example, the limit 05 — > 0. 
We have r(a5) in the denominator but also the 
product r(o5-f 2;i-|-Z4)r(— zi)r(— Z4) which is sin- 
gular in this limit. To reveal the singularity we 
take residues at zi = and Z4 = and obtain a 
factor r(o5) so that the limit becomes nontrivial. 

In the second step, one resolves the singularity 
structure in e, taking residues and shifting con- 
tours, with the goal to obtain a sum of integrals 
where one can expand integrands in Laurent se- 
ries in e. One can apply two strategies formu- 
lated in ^ and [2]. According to the first strat- 
egy, one performs an analysis of the integrand to 
reveal how poles in e arise. The guiding princi- 
ple is that the product r(o + z)T{h — z), where o 
and h can depend on the rest of the integration 
variables, generates, due to the integration over 
z, the singularity of the type r(a -I- h). So, one 
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thinks of integrations in various orders and then 
identifies some 'key' gamma functions which are 
crucial for the generation of poles in e. Then one 
takes residues and shifts contours, starting from 
first poles of these key gamma functions. For con- 
tributions of the residues, the same analysis and 
procedure is appfied. (See [2^1 for details.) 

Within the second strategy 0, one chooses an 
initial value of e and values of the real parts of the 
integration variables, Zi,w, . . . in such a way that 
one can integrate over straight lines. Then one 
tends e to zero and whenever the real part of the 
argument of some gamma function vanishes one 
crosses this pole and adds a corresponding residue 
which has one integration less and is treated as 
the initial integral within the same procedure. 

The third step of the method is to evaluate in- 
tegrals expanded in e after the second step. Here 
one can use Barnes lemmas and their corollaries 
to perform some of the MB integrations explic- 
itly. In the last integrations which usually carry 
dependence on the masses and kinematic invari- 
ants, one closes contour in the complex plane and 
sums up corresponding series. 

3. Applying Grobner bases to solve IBP 
relations 

Let A be the ring of polynomials of n variables 
xi, . . . ,Xn and Z C ^ be an ideal with a basis 
{/ij /2, ■ • ■ , /fc}^- A basis is a Grobner basis if 
any polynomial g € X is reduced modulo this ba- 
sis to zero for any sequence of reductions. To 
define reduction one needs an ordering of mono- 
mials cx]^ . . . x^" . For the lexicographical order- 
ing, (ii, . . . ,i„) >- (ji, . . . , j„) (first monomial is 
higher than the second one) if there is Z < n such 
that ii = ji, 12 = j2, ■ • • , ii-i = ji-i and ii > ji. 
Then the leading term g of a polynomial g is the 
monomial which is higher than any other mono- 
mial. Now, the reduction is defined as follows. 
Suppose that the leading term of a given g is divis- 
ible by the leading term or some /i, i.e. g — Qfi- 

ring is a set with multiplication and addition. A subset 
T of a ring R is called a left (right) ideal if (i) for any 
a,b £ 2 one has a + b G I and (ii) for any a £ X, c g /? 
one has ca £ X (ac € X respectively). A subset {fi} of X 
is a basis if any g (z X equals ^ r^/i for some € A. 



Let gi = g — Qfi- The leading term of gi is lower 
than the leading term of g and gi G X if and only 
if g G X. One can continue and proceed with 
gi as with g and obtain similarly g2, gs, ■ ■ ■■ The 
procedure is repeated until one obtains gi = or 
an element gi such that gi is not divisible by any 
leading term fi. One says that g is reduced to gi 
modulo the basis {fi, f2, ■ ■ ■ , fk}- 

If a given basis is not a GB one can construct 
a GB starting from it and using the Buchberger 
algorithm. Suppose that fi — wqi and fj — wqj 
where w, qi and qj are monomials and w is not a 
constant. Define S'(/i, fj) — fiqj — fjqi ■ Reduce it 
modulo {fi} as described above. If the reduction 
gives a non-zero polynomial add it to the initial 
basis as fk+i- Consider then such S'-polynomials 
for other pairs of elements (including the new el- 
ement) and reduce them modulo the 'current' ba- 
sis. If there is nothing to do one obtains a GB. 
It has been proven by Buchberger 9 that such 
procedure stops after a finite number of steps. 

A classical problem is to find out whether a 
given element g € A is a member of X or not. 
This problem can be solved by choosing an order- 
ing and constructing the corresponding GB with 
the help of the Buchberger algorithm. After that, 
one applies the reduction procedure modulo the 
constructed GB to verify whether a given element 
belongs to the given ideal X. If g is a monomial 
and {51, 32, ■••,<?;} is a GB, then 

Xi^ . . . a;„" = ^ ^ rigi + ^ ^ Cij^^,,,^i^Xi . . . x^ , (2) 

1=1 

where none of the monomials in the second 
sum is divisible by highest terms of gi. Since 
{gi, 32, • ■ • , 5;} is a GB, the last sum vanishes if 
and only if g g X. If g does not belong to I, the 
last sum is non-zero, and the part of g belonging 
to X is completely included in the first sum. 

Let us now denote by X the (left) ideal gen- 
erated by the elements fi which define IBP re- 
lations. To solve IBP relations is to express the 
value of F at an arbitrary point (ai, 02, . . . , a„) 
in terms of the values of in a few spe- 
cially chosen points, i.e. master integrals. This 
problem can be solved similarly to the alge- 
braic problem described above. Let us think 
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of the case, where ah the indices are posi- 
tive. Theni^(ai,a2,...,a„) = {Y^'-\ . .Y^:--^ ■ 
1, . . . , 1) . In this case it is reasonable to con- 
sider the operators Yi as the main operators and 
get rid of the operators Y~^ by multiplying (of 
course, from the left) the operators fi by suffi- 
ciently large powers of the operators Yi. Let us 
assume that we are dealing with such fi . 

Let us observe that the situation is quite similar 
to the above algebraic problem: instead of poly- 
nomials in the variables xi, . . . ,Xm we have poly- 
nomials in the shift operators Yi , . . . , y„ . The 
natural idea is to turn, from the initial basis 
{/i, ■•■,//}, to a GB {gi,...,gv}. Indeed, it 
is known that this can be done similarly to the 
above case: one introduces an ordering and the 
notion of the highest term which define the re- 
duction modulo basis, then one can apply a gen- 
eralization of the Buchberger algorithm. The mo- 
tivation is the same: this is the GB that charac- 
terizes the given ideal in the 'best' way, so that 
the parts which belong to the given ideal do not 
belong to the second sum in Eventually, 
one obtains a similar relation, Y^^^^ . . . Y^"^"^^ = 

E!=i rigi + Y.Czi,...,z„Y^'~^ ■ ■ ■ Y,]"-'^ , Let us ap- 
ply it to -F, take the value at = 1 and use the 
fact that the operators of I give zero on F. We 
obtain F(ai,...,a„) = Y.Cit.,...,ir,F{ii,...,in)- 
Integrals on the right-hand sides of such relations 
(for various ai, . . . , a„) are master integrals. 

However, any implementation of this similarity, 
using a generalization of the classical Buchberger 
algorithm, meets a lot of difficulties. The given 
problem of solving IBP relations becomes much 
more complicated at least because one has to con- 
sider also non-positive indices a^. 

Another complication is the presence of 
the variables ai as non-commutative operators. 
Moreover, coefficients in monomials in the shift 
operators, Yi^ can vanish at some points. Still let 
us imagine a situation where one can apply the 
Buchberger algorithm to construct a generaliza- 
tion of the GB for solving a reduction problem 
in the case of positive indices. Simplest examples 
show that the number of the master integrals as- 
sociated with this region can be greater than the 
number of the 'true' master integrals. 



These complications lead to the natural idea 
[T^ to change the strategy based on GB. For a 
given family of integrals, F{ai, . . . , a„), the whole 
region for each is decomposed into the regions 
with ai > and < 0. The whole region of the 
multi-indices is decomposed into 2" sectors a^, = 
{(fli, . . . , a„) : fli > if i G I' , < if « ^ ;/} 
labelled by subsets C {1, . . . , n}. In a given 
Oij it is natural to consider Yi for i ^ v and Y^^ 
for other i as basic operators. Another important 
point in extending Buchberger algorithm, is tak- 
ing into account boundary conditions, i.e. specify 
the (trivial) sectors where the integrals vanish. 

One has to construct a basis of Grobner 
type for each non-trivial sector cr^,. The basic 
operations are the same as above, i.e. calculating 
S'-polynomials and reducing them modulo current 
basis, with a chosen ordering. The goal is to con- 
struct a so-called sector basis IT (s-basis) which 
provides the possibility of a reduction to master 
integrals and integrals of lower sectors, i.e. a,yi 
for v' C ly. This point of the strategy is based 
on multiple examples of solving IBP relations by 
hand, where one tried to reduce indices to zero. 
It turns out that, within this strategy, one would 
construct a true GB only in the case of the sector 
(70 but this sector is always trivial. 

After constructing s-bases for all non-trivial 
sectors one obtains a recursive (with respect to 
the sectors) procedure to evaluate F{ai, . . . , a„) 
at any point and thereby reduce a given integral 
to master integrals. (See for details of the 
algorithm.) Examples have shown that the algo- 
rithm works at the level of modern calculations. 
In jl2| . a non-trivial example of integrals with 
seven indices was considered and, in [53], reduc- 
tion problems for two families of HQET integrals 
with nine indices were solved successfully. 

4. Some perspectives 

For families of complicated Feynman integrals, 
a reduction procedure was considered obligatory. 
On the other hand, it was clear that at least 
the second strategy of resolving singularities in 
e within multiple MB representations could be 
formulated algorithmically and implemented on 
a computer. (I believe that the first strategy can 
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also be automated and that the two strategies can 
be combined in order to achieve an optimization 
of calculations.) Recently, two algorithmic for- 
mulations have appeared j26l27| . so that, in com- 
plicated situations where one fails to solve the 
reduction problem by Laporta's algorithm, or by 
Baikov's method, or by using GB, one can try 
to calculate every integral, at least numerically, 
using these formulations. 

In fact, the second of these algorithms TT has 
been already implemented in Mathematica. This 
code provides a very good precision. For the ten- 
nis court integral discussed above this numerical 
integration provides excellent agreement with the 
analytic result of 17 . (If such algorithm existed 
a year ago, the authors of 17 would be satisfied 
by this powerful check and would not calculate 
asymptotic behavior when s/t — > by expansion 
by regions |28| :-) ) Anyway, evaluating compli- 
cated Feynman integrals, without reduction, be- 
comes now a reliable alternative. 
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